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Abstract 

We consider the Lawrence-Doniach model for layered superconductors, in which 
stacks of parallel superconducting planes are coupled via the Josephson effect. We 
assume that the superconductor is placed in an external magnetic field oriented parallel 
to the superconducting planes and study periodic lattice configurations in the limit 
as the Josephson coupling parameter r — > 0. This limit leads to the "transparent 
state" discussed in the physics literature, which is observed in very anisotropic high-T c 
superconductors at sufficiently high applied fields and below a critical temperature. 
We use a Lyapunov-Schmidt reduction to prove that energy minimization uniquely 
determines the geometry of the optimal vortex lattice: a period-2 (in the layers) array 
proposed by Bulaevskii & Clem. Finally, we discuss the apparent conflict with previous 
results for finite- width samples, in which the minimizer in the small coupling regime 
takes the form of "vortex planes" (introduced by Theodorakis and Kuplevakhsky.) 
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1 Introduction 

In this paper we continue the analysis of layered superconductors in a parallel external 
magnetic field started in our previous paper [[AlBcBr 00|| . In the previous paper we considered 



the case of a superconducting sample of finite width; here we treat periodic solutions in an 
infinitely wide superconductor. 

The Lawrence-Doniach model ||LaDo Tl]] is a mesoscopic Ginzburg-Landau type model 



for superconducting materials with a planar layered structure, and was originally applied 
to study organic superconductors and other superconducting composites manufactured by 
deposing successive thin layers of superconducting metal with interposing insulating films. 
Interest in this model has been spectacularly revived by the discovery of high temperature 
superconductors, since nearly all of these materials exhibit a distinctly layered structure. 
Indeed, a pure mono crystalline sample of cuprate high-T c material (such as Bi 2 Sr 2 CaCu 2 0s 
(BSCCO), Tl 2 Ba 2 CaCu 2 8 (TBCCO), or to a lesser extent YBa 2 Cu 3 7 (YBCO)) consists 
of copper oxide superconducting planes stacked with intervening insulating (or weakly su- 
perconducting) planes. 

The Lawrence-Doniach model. We model the layered superconductor as an infinite stack 
of superconducting planes, each parallel to the xy-pl&ne and with uniform separation distance 
p. The planes are thus described by those (x, y, z) with (x, y) e R 2 and z = z n := np, n E Z. 
We impose an applied "external" magnetic field, of constant magnitude H and lying parallel 
to the planes, along the y-direction, H = Hy. We make the ansatz that the local magnetic 
field inside the sample will be everywhere independent of y and point in the y-direction, 

h(x,y,z) = h(x,z) y. 

The vector potential A may then be chosen to lie in the xz-plane, 

_ ( dA dA \ 

A(x, y, z) = A x (x, z) x + A z (x, z)z, h = curl A = ( - J y. 

We study configurations which are bi-periodic, in the sense that the physically observable 
quantities (the magnetic field, the density of superconducting electrons, and their currents) 
are doubly periodic functions in the xz-plane. More precisely, given q > 0, iV e N, and 
s G R we define the fundamental domain of periodicity to be the parallelogram II = TlN,s,q 
spanned by the vectors t\ = (2q, 0) and e 2 = (s, Np): 

II := {(x, z) = t x e x + t 2 e 2 , < t u t 2 < 1}. 



In other words, our fields and currents will be 2g-periodic in x (along each superconducting 
plane), and will repeat themselves after each N planes but with a horizontal translation of 
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s. We denote 



np and x n = ns/N, n G Z, 



so the part of the nth SC plane which lies within the basic period module II is described by 
the points (x, z n ) with x n < x < x n + 2q. 

In each plane we define a (complex- valued) superconducting order parameter ip n (x), 
n G Z. We choose units in which \ip n \ = 1 represents a purely superconducting state. With 
the assumption that the physically observable quantities are H-periodic we may measure the 
free energy over a single period to obtain the following functional, 




x n +2q 



r N rx n +2q 

+ 2 P ^ J 

4 . 1 Jx n 



— iA x ip n 

dx I 



1p n - V'n-iexp i 





AJx, s) ds 



dx (1) 



dx 



H dxdz 



Note that is expressed in non-dimensional units, chosen such that the in-plane pene- 
tration depth \ ab — 1, k = Xab/iab is the Ginzburg Landau parameter, r is the interlayer 
coupling parameter (or Josephson coupling parameter,) and the magnetic fields are measured 
in units of H c /k, where H c is the thermodynamic critical field. (See [ [I 1 96|| .) 

The coupling between the superconducting planes given by the second sum in Q^ er sim- 
ulates the Josephson effect, by which superconducting electrons travel from one supercon- 
ducting region to another by quantum mechanical tunnelling. We will see this explicitly in 
the Euler-Lagrange equations, where the currents in the gaps between planes will be deter- 
mined by the sine of the gauge-invariant phase difference. The interlayer coupling parameter 
r gives the strength of the Josephson coupling. 

As in our previous study of the boundary-value problem for the Lawrence-Doniach sys- 
tem fAlBeBr "00|1 we will study the minimizers (and low-energy solutions) for r near zero. 
Indeed, for the high-T c cuprates at a temperature sufficiently below T c we expect r is a small 
parameter in the Lawrence-Doniach model. In particular, for BSCCO or TBCCO we expect 
r(0) ~ 1CT 4 -10~ 3 . (See "Physical background" below.) 

Results. We find that for r ~ there is an unique periodic solution of the Lawrence- 
Doniach system which attains the minimum energy (per unit area) among any other periodic 
configuration, with any period geometry. This absolute minimizing solution has fundamental 
domain ILv iSj(J with N — 1, s — q — n/Hp: in other words, the currents and field are 2tt /Hp- 
periodic in x, and shifting from one plane to the next results only in a horizontal translation 



S. Alama, A.J. Berlinsky, & L. Bronsard 



4 



by a half-period s = 7c/Hp. For example, the magnetic field satisfies 

h(x + z \ = h(x, z) = h ( x + ——, z + p ] . (2) 



Hp 1 

Clearly this condition implies that the magnetic field, superconducting electron density and 
currents are 2p-periodic in z. In addition, the flux per period for this solution is exactly 
one quantum fluxoid, 2tc. This configuration is optimal in the following sense: each choice 
of period geometry and quantized flux determines a function space in which to minimize 

'• The configuration described above attains the minimum of free energy per unit cross- 
sectional area among all possible choices of U-N, s ,q and flux quantization. In fact our result 
says more: for any given geometry and flux, then for all r sufficiently small the minimum 
energy per unit cross-sectional area in that class is either always strictly larger than the 
_2_-bi-periodic solution above, or exactly equal. The energies coincide if and only if 
the given lattice ILv s is commensurate with IL _s_ and its minimizer coincides with the 
(unique) III, -minimizer described above. (See Theorem |4.1| .) 

The minimization among configurations and among geometries and fluxes is entirely the 
result of direct rigorous analysis of the Lawrence-Doniach system, without recourse to nu- 
merical approximation or further restrictions other than the periodic ansatz. This is not 
like the case for the Abrikosov lattice in the Ginzburg-Landau model near H c2 , where the 
energy-minimizing geometry and flux quantization were determined numerically by com- 
parisons within a finite collection of configurations. Our result confirms the prediction of 
Bulaevskh & Clem [|BuCm iJI ], who claimed that for sufficiently large applied fields H the 
energy minimizers should form a period-two vortex lattice satisfying an ansatz of the form 

©■ 

We also consider the case of finitely many superconducting planes, with periodic observ- 
able quantities. The result we obtain is much the same as in the bi-periodic case: mini- 
mization of the free energy per period strip dictates the choice of period 2q = j^, and the 
minimizing solutions are (to order r) 2p-periodic in z in the interior of the sample. (See 
Theorem |7.1| . ) 

In both settings we obtain a "transparent state" as observed in experiments on BSCCO 
by Kes, Aarts, Vinokur, & ver der Beek |[KAVB 90| (see "Physical background" below.) 
Superconductivity is essentially unaffected by the presence of strong magnetic fields, \ip n \ = 
1 — 0(r), while the magnetic field is virtually unscreened by the presence of the supercon- 
ducting planes, h = H + 0(r). In particular, the order parameters are never zero, and the 
"vortices" , which correspond to the local maxima of the local magnetic field, fit entirely in 
the gaps between the superconducting layers. In the physics literature these are referred to 
as Josephson vortices as opposed to the Abrikosov vortices characterized by the vanishing 
order parameter in the Ginzburg-Landau model. 
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The juxtaposition of the results obtained for the finite- width samples in our previous 
paper ||AlBcBr 00|1 with the result stated above may seem contradictory: in |AlBeBr OOH we 
prove that for any sample of finite width — L < x < L (and for r ~ 0) the unique absolute 
energy minimizer occurs when the Josephson vortices are aligned vertically in a "vortex 
planes" geometry. These solutions were described by the physicists Theorodakis |[Th 90|1 and 
Kuplevakhsky ||K 99|| , the latter claiming that they were the only solutions to the Lawrence- 
Doniach system. Our explanation for this duality rests on the dependence of the interval of 
validity of the perturbative regime r ~ on the length of the sample size: the larger the 
sample, the smaller the interval in r for which the results of ||AlBcBr 00|| are valid. We will 
provide a more complete comparison of the two solutions in section 8. 



Physical background. A first attempt to model layered superconductors superconductors 
is by an anisotropic Ginzburg-Landau model, which treats the sample as a three-dimensional 
solid with anisotropic material parameters. In particular, the coherence length along the 
perpendicular to the planes £ c is assumed to be different from the value of the coherence 
length £ a b within the planes (see [ T 96| | or |CDG 95 ], for example) and the anisotropy is 
measured by an "effective mass ratio" T := ^&/Cc- For example lye | |1 92| | cites values of 
Fybco — 49, T B scco — 3025, and T TBC co > 90000. For certain materials (such as YBCO) 
and temperatures close to the critical temperature T c this approximation seems valid, but 
for more anisotropic superconductors the anisotropic Ginzburg-Landau model does not give 
a good qualitative or quantitative description of experimental observations. 

Experimental evidence of the failure of the anisotropic Ginzburg-Landau model was 
observed by Kes, Aarts, Vinokur, and van der Beek | KAVB 90 1, in the situation where 
the sample is placed in a strong magnetic field oriented parallel to the superconducting 
planes. Experiments reveal a crossover between "three-dimensional" behavior (governed by 
the anisotropic Ginzburg-Landau model) and "two-dimensional" behavior when the tem- 
perature is lowered beyond a critical value T c0 < T c . The "two-dimensional" regime is 
characterized by a "magnetically transparent" state, in which the applied magnetic field pen- 
etrates completely between the planes, virtually unscreened by the superconductor. These 
observations are not in agreement with the anisotropic Ginzburg-Landau model, where the 
magnetic field is largely expelled from the bulk and isolated vortices appear in a triangu- 
lar "Abrikosov lattice." For highly anisotropic materials (such as BSCCO and TBCCO) 
the two-dimensional regime occurs within one degree Kelvin of T c , and an analysis of the 
transparent state requires a model which addresses the discrete nature of the material: the 
Lawrence-Doniach model. 

To motivate the treatment of r as a small parameter, we re-write r in the original dimen- 
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sional coordinates, 

r = r{T) 



XjP 2 XjK 2 p 2 ' 

where A j = Xj \ a b, p = p \ ab are the actual physical Josephson penetration depth and 
separation distance. The length scale Aj in the Lawrence-Doniach model is directly related 
to the effective mass ratio of the anisotropic Ginzburg-Landau model via T = Aj. The 
temperature dependence of r(T) is determined by the conventional dependences of ~ 
Co&Wfrfy anc ^ ^-ab — ^ob(O) t T -t " ^ or ^ ne hiS n ~^c cuprates, £ a ft(0) and p are of the same order 
of magnitude, and hence for highly anisotropic superconductors at temperatures sufficiently 
below T c , r ~ A J 2 = T" 1 is small. In particular, for BSCCO or TBCCO we expect r(0) ~ 
10- 4 ,10" 3 . 

We note that Chapman, Du, & Gunzburger ||CDG 95| 1 have proven that solutions of the 



Lawrence-Doniach model converge to solutions of the anisotropic Ginzburg-Landau model 
(and in particular the convergence of energy minimizers) under the limit p — > with k, T 
fixed. This limit does not correspond to our "two-dimensional" regime, since it would send 
r — > oo and is therefore far from the "weak coupling" of superconducting planes observed 
KAVB 90 . Indeed, condition (1) in | ]KAVB 9l| , which defines the dimensional cross-over 



in 



point T = T c q, is equivalent to r(T c o) < 1 in the Lawrence-Doniach model. On the other 
hand, by fixing \ ab = 1 in our units the limit T — > T c (from below) sends p = p/\ a b(T) — > 0. 
Therefore it is not surprising to recover "three-dimensional" behavior for p — > 0. 

Methods. As in our previous study of the boundary-value problem for the Lawrence- 
Doniach system ||AlBeBr 00[| we will use a degenerate perturbation approach to study the 



minimizers (and low-energy solutions) r near zero. However, the periodic and bi-periodic 
settings are more subtle and present some additional complications in applying the method. 

First, periodic magnetic fields and currents are generally represented by non-periodic 
order parameters and potentials (<p n ,A). In order to define a variational setting for the bi- 
periodic problem we adopt the idea of 't Hooft ||tH 79fl to define spaces of functions which are 



periodic up to gauge transformation. (See section 2.) These spaces are rather complicated 
from the functional analytic point of view, and therefore it is important to find an equivalent 
formulation. We do this by making a gauge change (in section 3) to arrive at a family of 
affine Hilbert spaces representing the periodic configurations. 

For the finite- width case considered in | |AlBcBr 0(J| j the crucial observation was that when 



r = the planes decoupled, and the energy could be minimized explicitly. Even after gauge 
symmetries were removed the r = problem exhibited an additional symmetry, correspond- 
ing to an (A^ — l)-dimensional torus action (where denotes the number of superconducting 
planes,) and the minimization problem at r = degenerated on a finite dimensional hyper- 
plane in function space. When r ^ this symmetry was broken and by a Lyapunov-Schmidt 
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decomposition we reduced the problem of finding solutions with r ~ to a finite dimensional 
variational problem on this hyperplane. In this finite-width setting, the minimum value of 
energy was 0(r). 

In the bi-periodic case the r = problem already dictates a choice of period and flux: 
the minimum energy for the configuration , s>q will be 0(r) if and only if q = j^, m G 
N, and the flux through IIv.s,? is exactly 2TrmN. (See Lemma [4.2| .) Within this more 
restricted class of geometries the r = problem can again be solved explicitly, and the 
same degenerate perturbation method as in fAlBeBr "00[| applies to calculate the minimum 



energy as an expansion in r for < r << 1. Unlike the finite- width case, this expansion 
yields a constant term at order r, and the computation must be carried out to order r 2 in 
order to select a point on the degenerate hyperplane which minimizes energy with r ~ 0. 
The resulting finite dimensional minimization problem may be completely resolved to give a 
clear choice among the remaining parameters. (See Theorem |6.1| .) 

The case of finitely many layers is similar (and simpler) than the bi-periodic case; we 
present a sketch in section 7 of the modifications necessary to treat that setting. 

As in ||AlBeBr "DO] we recognize that in a real superconductor the parameter r, while small, 



is not infinitessimal, which raises the question of the validity of asymptotic results in the 
regime r ~ for physically relevent parameter choices. In section 5 of ||AlBeBr 00|| we derived 



a lower bound on the radius of validity of the r-expansion in terms of the other parameters 
in the problem. The same analysis applies in the periodic problem, and we conclude that the 
radius of validity is enhanced with smaller period q and k, and larger external field H. Since 
our optimal configuration ties the period q to the reciprocal of the field H, we conclude that 
the result is most applicable in strong external fields. In particular, the expansion of the 
periodic problem is more likely to be applicable than the expansion obtained in fAlBeBr "00f| , 
since the validity of finite-width expansions deteriorates for large samples, whereas treating 
the sample as infinitely large shifts the dependence to the period, which can be small with 
large H. 

We briefly outline the organization of the paper: in the second section we introduce a 
functional analytic framework for the doubly periodic problem, based on the elegant formu- 
lation of gauge periodicity due to ! t Hooft | |tH 73J . The third section defines an equivalent 



space, eliminating bothersome symmetries due to gauge and continuous translation invari- 
ance. In section 4 we show how minimization at r = forces the correct choice of period 
and flux, and the fifth section reviews the perturbation method used in [|AlBeBr 00|| . For the 



periodic problems it is necessary to expand the energy to order r 2 , and hence the solutions to 
the projected equations (coming from the Lyapunov-Schmidt reduction) must be calculated 
to order r: these computations occupy section 6. The seventh section sketches the procedure 
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for the case of finitely many planes, which is conceptually simpler. 

Acknowledgement. SA and LB wish to thank F. Helein and the Ecole Normale Superieure- 
Cachan for their kind hospitality during the Winter and Spring of 2000. 



2 Variational formulation. 



We begin with the doubly periodic case: we adopt the notation of the introduction, and 
assume that the gauge invariant quantities (fields and currents) are periodic in each SC 
plane, and allow the pattern of fields and currents could repeat themselves with a horizontal 
shift in x after each period in n. As before, for given N G N, s G R and q > we define the 
fundamental domain II = IIjv,s,g as in the introduction. We seek solutions to the LD system 
which are periodic with respect to the lattice generated by e*i, e 2 , 



h(x + 2q, z) = h(x, z) = h(x + s, z + Np), 



(3) 



and similarly for the other gauge-invariant quantities: the density of superconducting elec- 
trons \ip n \, the in-plane current density j x n \ and the Josephson current density in the gaps, 
j( n \ defined by 



J. 



ii n) := Im 



J 



(n) ._ rK P 



Im 



Of course when s = we would have simple periodicity in z, with period Np. 

The subtlety of the periodic problems is that periodic magnetic fields and currents are 
generally represented by non-periodic potentials A and order parameter if) n . One setting for 
such periodic problems is via t'Hooft boundary conditions [ tH 79|| , for which one demands 
that A and ip n be periodic up to a family of gauge transformations from one period cell to 
the next. Let N,q,s be fixed constants which determine the period lattice e*i, e*2 as above. 
We say that 

A)erR = W(N, s, q) (4) 

if ip n G Hi oc (R), A G Hl oc (R 2 ; R 2 ) and there exist two functions uj x and u z in if^R 2 ) such 
that: 



ip n (x + 2q) = ip n (x) exp[iu x (x, z n )], ip n+N (x + s) = ip n (x) exp[iu z (x, z n )], 



(5) 



A(x + 2q, z) = A(x, z) + Vw x (x, z), A(x + s, z + Np) = A(x, z) + Vu z (x, z). (6) 

Since ijj n+ ^{x-\-2q) can be related back to ip n (x) in two different ways (by commuting the 
order in which the two rules in ([|) are applied,) this imposes a constraint on the functions 
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u x ,u z which are permissible in this definition. For any configuration (ip n ,A) G tTi there 
exists an integer K e Z such that 

u x (x, z n ) - u x (x + s, z n+N ) + lu z (x + 2q, z n ) - u z (x, z n ) = 2nK. (7) 

Clearly, by continuity K may be chosen independent of x; we will see below that it is also 
independent of n. The important role played by K is revealed by integrating the magnetic 
field strength over the basic unit period cell, applying Stokes' Theorem, and substituting the 
identity (0): 

if cml A dxdz = A ■ ds (8) 

JJU JdU 

= u x (0, 0) - w x {s, Np) + uj z (2q, 0) - cu z (0, 0) 
= 2nK. 

Hence the total flux per period cell is quantized for any element of tTi.. Note also that 
the periodicity of h = curl A implies that the constant K in (]?]) may indeed be chosen 
independent of n. 

Since every (?/>, A) e tTi is associated to an integer K via a continuous map, 

(ip,A)^—JJ cmlAdxdz, 

this discrete choice separates the space tH into disjoint connected components, 

fH(N,s,q) = |J tH(N,s,q,K). 
Kez 

We may of course minimize the Lawrence-Doniach energy in each component tTC(N, s, q, K) 



separately, and indeed Lemma O below will indicate the optimal choice of K, for appropriate 
(N, 8 ,q). 



Remark 2.1 These periodic configurations may also be described within the context of 
complex line bundles over the 2-torus T, defined as R 2 modulo the discrete lattice generated 
by e*i, e*2. A configuration (ip n , A) defines a connection on T by interpreting A as a one-form, 
while the if> n (x) are restrictions to the iV planes in II of a section of a complex line bundle 
over T, with structure group U(l). The 't Hooft conditions incorporate the U(l) group 
action as we pass from one coordinate patch to another on the manifold T. In this setting, 
the magnetic field h = curl A is the curvature of the connection, and (f|) is a Gauss-Bonnet 
relation between the total curvature and the Euler number K, reflecting the fact that the 
bundle is nontrivial. 
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While this is a very elegant formulation of the periodic problem it poses some practical 
problems for analysis of the variational problem because of the introduction of the aux- 
iliary functions u x ,u z as well as the usual degeneracies associated with gauge symmetry 
Fortunately we will be able to simplify the setting of the problem in two ways. First, by 
fixing a gauge we find a much simpler Hilbert manifold setting. This will be done in sec- 
tion 3. Second, we will show (see Lemma [4.2| below), that the least-energy solutions will 
have |^> n (x)| ~ 1, and therefore it will be convenient to use polar coordinates for ip n in order 
to treat the phases more directly. In the remainder of this section we introduce subspaces of 
tH for ip n in polar form and describe their properties and the Euler-Lagrange equations for 
critical points in these spaces. 

To define the space tTL in terms of f n , 4> n , with ip n (x) = f n (x)e l ^ n ^ we must take 
into account that the periodicity conditions on <p n can only hold modulo 2it. For the second 
relation in (|5|) (translating in z from np — > (n + N)p) this is not important, since the constant 
factor of 2tt can be absorbed into u z . The degree of winding per period in the x-direction 
plays a particularly important role, and hence we assume that there exist integers k n such 
that 

<f) n (x + 2q) = (f) n (x)+uj x (x,z n ) + 2'Kk n , (j) n+N (x + s) = (j) n (x) + u z (x, z n ). (9) 

This definition does not uniquely determine the winding numbers k n , since adding the same 
constant multiple of 2tt to each merely reduces the value of the function uj x (x,z) by that 
same quantity. To remedy this problem we may (without loss of generality) fix the value 

k = 0. 

Returning to the calculation (0), we see that 

K = k n+N - k n = k N - k = k N , 

and hence the index K associated to the space tH(N, s, q, K) measures the net change in the 
winding number of n over one period Np in z. Of course the moduli of the order parameters 
is gauge-invariant, and is Il-periodic, 

f n (x + 2q) = f n (x) = f n+N (x + s). (10) 

We may now play the same game with the winding numbers k n (n 6 Z) in the polar 
representation (/„, n , A) as we did with the single index K for the configurations (ip n , A) e 
tTC. Since (/„, <f) n , A) determines (k n ) ne z via a continuous selection, the space of admissible 
configurations splits into disconnected classes. Fixing k := (fc ra ) n ez (with ko = and k n+ N — 
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k n = K) we define our basic space of bi-periodic configurations BV(N, s, q, k) by choosing 
the corresponding connected component of £H: 

BV(N, s, q, k) := {([f n , <f> n ),l) : f n , <f> n G H^R), A G ^ oc (R 2 ; R 2 ), and there 

exist uu x , lo z G H? oc {R 2 ) such that (03), ®, and © hold.}. 

Each component BV(N, s, q, k) describes a subclass of tTC(N, s, q, K) with K = fcjv, and 
therefore the flux quantization 

curl A dx dz = 2nK = 2-nk^ 

n 

holds for every configuration (f n ,<p n ,A) G BP. As for tTC(N,s,q,K), it will turn out that 
there are preferred classes BV(N, s, q, k) for which the free energy will be smallest possible. 
(See Lemma (O] below.) 

It is a simple calculation to verify that for any configuration (ip n , A) G tTC(N, s, q, K) 
(and for (f n ,<f) n ,A) G BV(N,s,q,k)) we recover the desired periodicity conditions (see ([|)) 
for the gauge-invariant quantities. Since only the gauge-invariant quantities appear in an 
expression of the free energy density, the energy density will be Il-periodic and hence we 
are justified in measuring the energy of the configuration over only one period cell. We may 
therefore define the Gibbs free energy Qf p for the configuration (ip n , A) G tH(N, s, q, K) as 
in the introduction. When ip n can be represented by polar coordinates f n , <p n we use the 
equivalent form: 



JV fX n +2q M 1 1 

n? P (/„, KA) = P £ / ~{f n - I)' + -(fn? + -(<f>n- A X {Xn, y)? f 

+ nP E / (fn + fLl ~ 2/n/n-l OOS^x)) dx 

1 n=l Jxn 



dx 



1 ff (dA x 8A Z r V , , 
fir JJu \ oz ox J 

where $ njn _i(:r) is the gauge invariant phase difference, 

®n,n-l{x) ■= <Pn{x) ~ <f> n -l(x) - / A z (x, z) dz. 

■> Zn-l 

We observe that $ ni „_i(x) is periodic modulo 2n, 

$n tn _i(a; + 2q) = $ njTl _i(a;) + 2n(k n - fc n _i). 
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Remark 2.2 Note that although we have defined (f n , <p n , A) G BV(N, s, q, k) globally in R 2 , 
in fact the gauge-invariant quantities f n , ® n>n -i, V n := (4>' n — A(x, z n ), and h(x, z) are globally 
determined by the values of f n (x), (p n (x), and A(x, z) for n = 1, . . . , N and (x, z) 6 II, and 
the values of u> z (x, 0). This is very clear for f n , V n , and h; for $ n>n _i, we require 4>q(x) to 
determine <$>i t o(x), and it is here that the 't Hooft condition comes into play via the function 
u z (x,0). In retrospect, we could have defined (f n ,(p n ,A) only on the period module II and 
uj x ,u z only in a neighborhood of the left and bottom edges of II. This would be enough to 
make sense of the free energy, and it allows interpretation of the 't Hooft conditions as true 
boundary conditions, 

ip n (x n + 2q) = i) n (x n ) exp[iu x (x n , z n )], n = l,...,N; 
^n{x + s) — ipo(x) exp[iu z (x, 0)], < x < 2q; 

Np 

A(x + 2q, z) = A(x, z) + Vu x (x, z), < x < s, z = x\ 

s 

A(x + s, Np) = A(x, 0) + Vu z (x, 0), < x < 2q. 



This is the interpretation of the 't Hooft conditions taken by Tarantello ||Ta 96|| for example. 



The Euler-Lagrange equations under 't Hooft conditions are virtually identical to those 
derived in [|AlBeBr 00|| except that we must take into account the periodic boundary condi- 



tions both in x and z. To obtain the equations we may choose smooth genuinely n-periodic 
test functions to vary each unknown in turn. We obtain: 

- -X + (fn ~ + - 2 ^'n ~ A X (X, Z n ff n (11) 

r 

= ~ [fn-l COS $ n>rl _i + f n+ i COS $„+!,„ - 2f n ) ] 

^(x,z) = z^z n , (12) 
h{x, z n +) - h(x, z n -) = -pfl{x){(j)' n - A x (x, z n ))\ (13) 
0^ = —^-fn{x)fn-i{x) sin$ ni „_i(x), for 2 n _i < z <z n . (14) 

These equations should be solved together with n7v iSjg -periodic boundary conditions for the 
observable quantities, f n (x), h(x,z), V n (x) := (4>' n {x) — A x (x,z n )), and sin $ nin _ 1 (a;). For 
later purposes, we also record the current conservation equation 

hlx _ A x( x > z n))) = ^(/n/n-isin$ nin _i - / n+ i/„sin$ n+1>n ). (15) 

This equation is not independent of the others: it can be derived by differentiating (pf) 
and substituting from ([LI]). This is not surprising since (fi n and A are related through gauge 
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invariance. We note that 

£ n) (x) = f n (x)(<j ) ' n (x)-A x (x,z n )) 
measures the current density within the n th superconducting plane, while 

2 

jt\ x ) = —^-fn{x)fn-l{x) Sin$ ni „_i(x) 

gives the Josephson current density in the gap between the (n — l) st and n th planes. In 
this way we may view fll5D as a semi-discrete version of the classical equation of continuity 
divj = 0. This is the conservation law corresponding to the £7(1) gauge invariance as 
guaranteed by Noether's Theorem. 

We also record a useful formula for which follows from Stokes' Theorem: 

$n,»-i(aO = l X (V n -V n -i)dx + p ( X h {n \x)dx + (0). (16) 



3 Fixing a gauge 

While the spaces tfH(N, s, q, K) and BV(N, s, q, k) are indeed unwieldy for the purposes of 
analysis, we will show that by fixing an appropriate gauge we can, without loss of generality, 
work in a much simpler setting. We begin by noting that the periodic problem has a larger 
symmetry group than the fixed interval problem studied in the previous paper fAlBeBr TKj . 
In addition to electromagnetic gauge invariance, 

fn^fn, (f) n (x) (f) n (x) - \(x,Z n ), A -> A - VA(x, z) , 

there is also translation invariance, 

(fn, 0n, A) -> (/„(■ - Xq), (j) n (- ~ Xq), A(- - X , •)), G R. 

While this last symmetry is only one-dimensional, for our degenerate perturbation theory it 
is convenient to eliminate all but the essential symmetries of the r = problem which are 
broken when r ^ 0, and so our choice of "gauge" will fix a translation as well. Note that 
there is also a discrete translation invariance in z, z — > z — kp, k G Z, but this symmetry 
(being discrete) will not create analytical difficulties for our method and hence plays a less 
important role. 

First consider the spaces BV(N, s, q, k) for which k ^ 0. (Recall that we have fixed 
ko = in the definition of BV(N,s,q,k).) By the discrete translation invariance in z we 
may relabel the z-axis if necessary in order to obtain 



h + = k . 



(17) 
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We observe that the average value of the local magnetic field is fixed by the choice of the 
space BV(N, s, q, k): indeed by @, 



is) 



where h = curl A. 

Next we define an appropriate space which eliminates gauge and translation symmetries. 
Assume again that k ^ 0. We say that (f n ,4> n ,A) G BV*(N, s,q, k) provided there are 
constants lj, d G R such that 



A G HUB 2 ; R 2 ), /„ G HUR), 4>n G #L( R 



(19) 



4> n (x + 2q) = (j) n (x) + uj + 2-nkn] 

k^TT 

<fi n +N{X + S) = <p n [x) + — — X + d] 

f n (x + 2q) = f n (x) = f n+N (x + s); 
l(x, z) = (h) (z, 0) + (d z £ , -d x , where 

f g //L(R 2 ) ? £(s + 2g, = + s, z + iVp) = z); J 

I 4> (x)dx = 0= / (j)i{x)dx. 

J—q J—q 

As usual, we may choose a definition of the phases so that 

< n (O) < 27r for each n G Z. 



(20) 
(21) 
(22) 

(23) 
(24) 

(25) 



In the case k = 0, we define BV*(N, s, q, 0) via (|i~9"|)-(p3|), and we demand only the vanishing 
of the mean value of <fto(x) in (with no restriction on the mean of 4>i(x)). 

It is easy to see that BV*(N, s,q, k) C BV(N,s,q,k). We must show that the two 
spaces are equivalent, in the sense that for any configuration in BV(N, s, q, k) there is a 
gauge transformation which sends it to an element of BV*(N,s,q,k). In particular, this 
implies that minimization of fl^ p in BV*(N,s,q,k) is equivalent to minimization of 
in BV(N, s, q, k). In addition, we will show that the choice of gauge in BV*(N, s, q, k) is a 
"Coulomb" gauge, in the sense that the H l (I\) norm of A is controlled by its curl in L 2 (II). 



Theorem 3.1 Let N, s, q, k be fixed. 
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(a) There exists a constant Co > such that for any (f n , <f) n , A) G BV*(N, s, q, k), 

\\A\\ 2 m{n) < Co||^||i,2(n), 

where II is the period module. 

(b) For any (/„, 0„, A) G BV(N } s, q, k) there exists A G Hf oc (H 2 ) and x G R stzc/i t/iat 
(/ n (x - x ), 0„(x - x Q ) - X(x - x , z„)> A(x - Xo, 2) - VA(x - x , z)) G BV*(N, s, g, fc). 

When = the choice of Xo in (b) is immaterial: we cannot remove translation invariance 

in that case. 

Proof of Theorem j| : 

Suppose (/ n ,0„,v4) G BV*, with A defined as in (|2"B|). Then A£ = h — (h) in R 2 , and 
therefore the estimate follows from standard elliptic regularity theory. 

To prove (b), assume (f n , (f> n , A) G BV(N, s, q, k) and let £ be a solution to the periodic 
problem 

A£ = curl A - (h), 

£(x + 2q, z) = £(x, z) = £(x + s, z + Np), 

and A = (<9 Z £, — d x £). By standard elliptic regularity theory £ G if 2 oc (R 2 ). Then curl (A — 
A) = so that there exists A G if 2 oc (R 2 ) with A = A — VA. In fact A is only determined up 
to a constant; set A = A — c, with c as yet undetermined. 
Define 

4>n ■= 4>n - X{x,Z n ). 

It follows that 

<p n (x + 2q) - n (x) = u x (x, z n ) - \(x + 2q, z n ) + X(x, z n ) + 2ixk n . (27) 
Since A is periodic in x, we have (for all z) 

Voo x (x, z) - VX(x + 2q, z) + VA(x, z) = 0, 
and hence there exists a constant of integration u so that 

u) x (x, z) — X(x + 2q, z) + X(x, z) = u> 
for all (x,z) G R 2 . Substituting this identity in fl2~T|) we verify ( pD|) for n . Similarly 
0„ +A r(x + s) - 0„(x) = ^ 2 (x, z n ) - A(x + S, 2„ + iV» + X(x, z n ), 

and since 

i(x + s, z + Np) - A(x, z) = (Np, 0)(h) 
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it follows that there is another constant of integration d such that 

u) x (x, z) — X(x + s, z + Np) + X(x, z) = Np(h) x + d 

for all (x,z) G R 2 . In particular using the quantization formula (|T8|) , it follows 

~ Kir 

<fi n +N{X + S) - <p n \x) = —X + d, 

and is satisfied. 

To complete the argument in case k 7^ , we define 

-t+q 



g{t) ■= J t (0i (a?) - <t>a(x)) dx. 



Recalling ([H]), k = ^ k%, so (|20D implies that — 0o(a ; ) ~ ^j 21 ^ for \x\ large. Therefore 
— > ±oo as t — >■ ±oo, or vice- versa. Since g(t) is continuous there exists x G R so that 
g(xo) = 0. We now choose the constant term in A to be 

c=— (f) 1 (x)dx = — 4> (x)dXi 

Zq JXQ — q Zq JXQ—q 

and the conclusion of (b) follows with A the desired gauge change and xq the translation. 

In case k = 0, we cannot assert the existence of such a translation xq, but we can still 
choose c = y I- q (po( x ) dx to obtain (j) with zero mean in [—q, q}. 



The spaces BV*(N, s, q, k) are quite reasonable from the analytical point of view. First, 
note that each (fi n (x) is the superposition of a linear function (with slope u> + 2nk n ) with a 
2g-periodic function, and hence <t>' n (x) is 2g-periodic. Since 

1 rs+q 

to = 4>\{x\ + 2q) — (f)i(xi) — 2nki and d=— \ (j)^{x)dx, 

Zq 



-n 



by the conditions (pUP, ( pT| ) the quantity 

iV 

Ell- 



n=l 



controls the i?/ oc (R) norms of any finite collection of <fi n , n G Z. 

Indeed, for any choice of parameters BV*(N, s, q, k) forms an affine Hilbert space. The 
tangent space (of admissible variations to (/„, n , A) G BV^) is the space E = E(N, s, q) 
consisting of all (u n , v n , a) such that 
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(1) u n ,v n G Hl oc (R) (n G Z) and there exist constants lj, d G R such that 

u n (x + 2g) = u n (x) = u n+N (x + s), 
v n (x + 2g) = t> n (;r) + lu, v n+N (x + s) = v n (x) + d. 

(2) a = (d z £, -d a £) for some £ G Hf oc {R 2 ) with 

£(x + 2g, z) = £(ar, *) = £(z + s, z + Np). 

In view of the above remarks on controlling the extension of elements of BV*(N,s,q,k) 
beyond the basic period II, we may choose a norm on E of the form 

N 

\\(u n , v n , a)\\% = [||ttn||H-i([, T „„ T „4-2 ? ]) + lknllffi([,T„.^+2g]) + JJ [cut\ a] 2 dx dz. 

71=1 ^ 

It is clear that Vt pp is a smooth (C°°) functional on BV*(N, s, g, fc). Since any configu- 
ration in BV(N, s, q, k) can be associated to an element in BV*(N, s, q, k) (via gauge trans- 
formation and translation) with the same energy, the minimizer of Q pp in BV*(N,s,q,k) 
will also minimize Q pp in BV(N, s, q, k), and hence satisfy the Lawrence-Doniach system of 
equations with periodic conditions. 

Finally, we note that the same procedure may be used to fix a gauge in the space 
tfH(N,s,q,K) of doubly gauge-periodic configurations (if) n ,A), to obtain a nice space 
tTC*(N, s, q, K) which eliminates gauge invariance. Indeed, we say (ip n , A) G tH*(N, s, q, K) 
if: 

(a) V„ e Hl 0C (R) and A G #/ 0C (R 2 ; R 2 ); 

(b) A satisfies (f§); 

(c) there exist constants cj, d G R such that 



ip n (x + 2q) = tp n (x) exp[ia>], ^„ + jv(z + s) = Vn(z) exp 



i x + d 



A version of Lemma |Q| holds true for spaces tTC*(N, s, q, K): for any (ip n , A) G tTC(N, s, q, K) 
there exists A G -£^ 2 C (R 2 ) such that 

(Vn exp(iA(-, z n )),A- VA) g *H.(JV, s, g, K). 

And for any (■?/>„, A) G t?i*(N, s, q, K) the estimate of Lemma |3.1| (a) also holds. 



S. Alama, A.J. Berlinsky, & L. Bronsard 



18 



4 Minimization at r = 

We are ready to begin the process of solving for the lowest-energy periodic solution to the 
Lawrence-Doniach system. As we have already seen the periodic problems are indexed by 
several parameters: the number of layers in the period module N; the geometry of the 
period parallelogram IT represented by the horizontal period along the SC layers 2q and the 
horizontal shift s when advancing N planes in z; and the winding numbers k n , which we 

— * 

will denote collectively by k. It will turn out that the geometry of the vortex lattice will be 
completely determined by energy minimization! 

Denote the minimum free energy per unit cross-sectional area over all admissible lattice 

— * 

configurations in BV*(N, s, q, k) by 

e r (N, s, q, k) := inf j ^^"'^ = (/», tfn, A) E BV*(N, s, q, k) 
Our goal is to minimize twice: first for each fixed period geometry HN, s ,q and choice of 

— * — * 

winding numbers k, and then to find the parameters N, s, q, k which give the least energy 
per unit cross sectional area. The main result we will prove is that, for < r << 1, the 
periodic solution with the least free energy among all possible lattice geometries is a period- 
2p in z lattice, with period 2q = in x. We denote by Z the integer sequence, k n — n in 
the following statement of our theorem: 

Theorem 4.1 For any choice of (N, s, q, k) there exists f = r(N, s, q, k) such that either: 

-> IT 71 

e r (iV, s, q, k) > e r (l, — , — , Z) for all r < f; 

or 

e r (N, s, q, k) — e r I 1, — , — , Z I for all r <r, 

in which case the minimizers are achieved in BV*(N, s,q,k), and they coincide with the 
period-2p in z, period 2q — Jj^ in x minimizers in BV*(1, Z). 

To completely clarify the above statement, we note that the magnetic field and the su- 
percurrents \ associated to elements of the space BV*(l, Z) make a horizontal 

translation of a half-period q 1 :— when we increase z by its period p: 

f n (x + 2q x ) = f n (x) = f n+1 (x + q x ), V n (x + 2q x ) = V n (x) = V n+1 (x + q x ) 

h(x + 2q 1 ,z)=h{x,z)=h{x + qi,z + p), ) (28) 

$ n ,n-i(x + 2gi) - 2vr = $„,„_! (a;) = $ n +l,n( X + 9i)> 
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where we recall V n (x) = (</>' n (x) — A x (x,z n )). It is clear that these solutions are also 2p- 
periodic in z. These solutions have the period structure proposed by Bulaevskii & Clem 
PuCm 91] . In section 6 we will show that the associated local magnetic field h(x,z), in- 



plane and Josephson currents and density of superconducting electrons f n admit an 

expansion near r = of the form: 



h(x, z) = H + r (-l) n+1 |^ cos(Hpx) + 0(r 2 ), z n -i < z < z n , 

= r (-l) n+1 j^cos(Hpx) + 0(r 2 ), 

j( n \x) = \ r (— l) n K 2 psin(Hpx) + 0(r 2 ), 

fn(x) = 1 - | + 0(r 2 ). 

The proof of Theorem 4.1 will be accomplished in many steps, concluding at the end 
of section 6. The remainder of this section is devoted to the study of the r = prob- 
lem. Section 5 introduces a degenerate perturbation method based on a Lyapunov-Schmidt 
decomposition at r = 0, and section 6 solves the reduced problems arising from the decom- 
position. 



(29) 



is 



Choosing the period and winding numbers. The first step in proving Theorem [O 
to consider the case r = 0, in which case the superconducting planes decouple. We know 
from the previous treatment of the finite-width problem in fAlBeBr "00] that the case r = 



is analogous to the self-dual point k = l/y/2 in the Ginzburg-Landau model in the sense 
that at that point minimizers satisfy a first-order Bogomolnyi system, in addition to the 
(second-order) Euler-Lagrange equations. Indeed, we will see that the reduction to a first- 
order system of equations is only possible when the period q and the winding numbers k 
are chosen appropriately. For a generic choice of q, k the minimum of energy will be much 
larger. 

Lemma 4.2 For r = 0, 

(a) inf{C/jf- p ('0 Jl , A) : (i/) n ,A) G tH*(N, s,q, K)} = if and only if there exists m G N 
such that Hpq = itlqh and K = m^N . 

(b) Assume Hpq = m^-n and K = rrioN, itlq G N. For any r > we have 

inf{^ p (^, A) : (ip n , A) G tH*{N, s, q,K)}< ^Nqpr. 

Moreover, there exist constants r$, C > such that whenever Q^ p (ip n , A) < ^-Nqpr, 
and < r < r , then \ip n (x)\ > 1 — CV 1 / 2 > for all x, n. 
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(c) 



inf 

BV*(N,s,q,k) 



if and only if there exists m G N such that Hpq = m 7r and k n = m^n for each n G Z. 



We remark that the constants C, ro in (b) will depend on N, s, H, mo, K. 

Proof: We begin with (a), (c). Assume q = k n = m n, K = k N = m N, with m G N, 

and choose 



nO 



>n0 



a n + nHpx, ip. 



nO 



o i(f>n0 



A o {x,z) = (Hz,0), 



(30) 



for any a n G R. It is easy to see that (/, 
with g* p & n0 ,A ) = and n pp (f 



nO, H>nO, 



A ) G BV*(N, s, q, k), (Vta, A ) G m*(iV, s, q, K), 



Aq) = 0. 



To see the opposite implications, assume first that mi BV ^ N s q ^VL pp = 0. By a simple 
estimate (see Proposition 2.1 of [AlBeBr 00j| ) any minimizing sequence must have /„ — ► 1 
uniformly It is then easy to see that the energy controls the norm of the minimizing sequence 
in BV*, and the minimizing sequence converges strongly to (f n , <fi n , A) G BV*(N, s, q, k), and 
attains the minimum of zero energy. Being a sum of positive terms, fi pp (f n o, 4>no, A ) = 
implies that each term is individually zero, 

f n = 0, <j)' n - A x (x, z n ) = 0, curl A — H. 

— * — * 

The vector potential A = (Hz, 0) is uniquely determined by the gauge given by BV*(N, s, q, k), 
and integration of the second relation leads to the solution space described by ( |3~Q|) above, 
with a n G R undetermined constants of integration. Therefore the question reduces to de- 
termine when the family fl3"0|) belongs to BV*(N, s, q, k). Applying condition (|20| ) we obtain 

2Hpqn — uj = 2nk n , for all n G Z. 

Since we have defined ko = 0, it follows that uj = and therefore Hpq = mir and k n = mn 
for some m G Z. The exact same argument applies to prove the analogous statement in (a). 

To prove (b), first observe that the upper bound on the infimum follows from inserting 
the test configuration ip n , A above into Q pp . Then we use the simple inequalities, 



d 

dx 
to obtain 

N 



(x) < 



' d_ 

dx 



iA x (x,z n ) ip r , 



x 



2\2 



;i-i^)ir<(i-iwr) 



N N r x n +2q \ 1 

p£l|iHvgill M w.+2 9 ] < « 2 pE/ ( 1 -w 2 ) 2 + - 5 

n=l n=l x " ™ 



ax 



47T -» 

< — 2 Q PP (^A)<1Nqpr. 
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By the Sobolev embedding theorem in one dimension we conclude that for r small enough 
there exists a constant C with |-0 n (a:)| > 1 — CV 1//2 and (b) is proven. 



Lemma ^72| justifies certain restrictions on the family of spaces BV*(N,s,q,k) when 
seeking the absolute minimizers of energy. In particular, in the remainder of the proof we 
fix 

7T -> 

Q — Qm '■— m ^ m G N, k = mZ. 

Hp 

Furthermore, statement (b) justifies the treatment of the problem in terms of the polar coor- 
dinates i/) n (x) = f n (x) exp(i(f) n (x)), which in turn allows us to restrict the winding numbers 
k n as above. In the space BV*(N, s, q m , mZ) the applied flux per period per plane is 27im, 
and from (||) the total flux through a period cell is 

curl A dzdx = 2nNm. 

n 

Therefore, for any configuration in BV*(N,s,q m ,mZ) the mean magnetic field in a period 
cell coincides with the applied magnetic field: 

C 1 ) = n — at~ 11 curl A dzdx = H. 



2q m Np JJn 

This suggests that we are indeed in the "transparent" state referred to in the Introduction. 

Note in addition that the choice m = 1 gives the smallest q for which the minimum energy 
is 0(r) for r small, and 2q% corresponds to the minimal period of the explicit minimizer (^0|) 
of the r = problem. Furthermore, the condition k = Z will imply that there is one 
Josephson vortex in between each adjacent pair of SC planes, per period in x. We will see 
that the minimizing configuration in BV*(N, s, q m , mZ) will in fact have minimal period 2q\, 
and hence reside in the space BV*(N, s, qi,Z). 

In summary, with the above restrictions on the period and winding numbers, the spaces 
of interest are 

/ vrm \ 
* \ N,S, Jlp~ ,mZ ) ' 



We recall that this space is defined as all (/ n , <f> n , A) satisfying (|T^), (p2|), (|2"3|), and ( p^) with 
q = q m and (h) = H, and 

4> n (x + 2q m ) = 4> n (x) + 2nmn; (31) 

(p n+N (x + s) = (j) n {x) + NpH x + d. (32) 



Identifying the solution set at r = 0. Having restricted our choice of period and winding 
numbers, thanks to Lemma |4]^ we may now identify the manifold of all minimizers to the 
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r = problem, and verify that the setting is appropriate to apply the degenerate perturbation 
theory as developed in [|AmCzE 8711 , fjAmBa 98| , [|AlBeBr 00| . 



Proposition 4.3 Assume q = q m , k = m G Z for some m G N, and r = 0. 

(a) inf {QQ P (f n , <p n , A) : (f n ,(j) n ,A) G BV*(N, s, q m , mZ)} = 0. The minimum value is 
attained, and the set of all minimizers is the (N — 1)- dimensional hyperplane 

S '■= {(fnAn,A)eBV*(N,s,q m ,mZ): f n = 1, 4> n (x) = a n + nHpx, (33) 
A x = Hz, A z = 0, where cto = = a%, (0:2, • • • , «jv) G R w_1 .} 

(b) For any element o = (/°,0°,A ) G S, the linearized operator D 2 Qq P (a) : E — > E 
defines a Fredholm operator with index zero. Moreover, 

T a S = kerD 2 n pp (a) ~ R^" 1 . (34) 

Recall that the space E = E(N, s, is the tangent space to BV*(N, s, q m , mZ), and was 
defined at the end of section 3. We define the second variation of energy as a quadratic form 
on E: in particular, for (/°, <f>° n , A ) G S and {u n , v n , a) G E, 

D 2 n o (/°),0°,^ O )K,^,a] = pJ2 \2ul + -[u' n } 2 + -[v' n -a x {x,z n )] 2 \dx 

n= l Jx n -q m IKK ) 

+ ff Icurlal 2 dxdz. 



k? JJn 



The proof of Proposition |4]^ (a) follows easily from the proof of Lemma |4.2j (c), and part 
(b) is clear from the form of D 2 Qo(f®), 0°, A ) above (see Proposition 2.1 in [[AlBeBr 00|| for 
details.) 

Note that T = T a S is independent of a 6 5, and that we may treat S as a (compact) 
(N — l)-torus T^ -1 , since the energy is 27r-periodic in each a n . Applying the condition (f2~T|) 
to an element of S we obtain the periodicity conditions for a n and 5 n := a n — a n -\. 

ajv+n — a n = d — Hps{n + N), 5 n+ N — S n = —Hps (mod 2ir), (35) 

where d is the constant which appears in the 't Hooft-type periodicity condition ([32|). Since 
we are given a = = a\ in the definition of the space BV*(N, s, q m , mZ), ( p5|) provides the 
relation aw = d — HNps. This reveals how a choice of («2, . . . , oln) determines all a n : 
determines the constant d, and the first equation of ( p5|) then generates all the others in the 
sequence. 

We also note that the a n (n = 2, . . . , N) may be used as parameters for the manifold S. 
However, it will be more convenient in the end to parametrize S by the phase differences, 

6 n . cx n OL n —\. 
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Since the values of (62, ■ ■ ■ , Sn) determine (a 2 , . . . , a at) uniquely, this is an equivalent choice 
of parametrization. We abuse notation and denote 

cr = a(5 2 , . . . ,5 N ) E S, 



5 Degenerate perturbation theory. 

We now perturb away from the degenerate minima of using a variational Lyapunov- 
Schmidt procedure, just as in [[AlBeBr 00[1 . This method has been used by Ambrosetti, Coti 
Zelati, & Ekeland ||Am(JzE 87|| , Abrosetti & Badiale ||AmBa 98|| , Li & Nirenberg |[LJN 98 



(and many others) in a variety of situations involving heteroclinic solutions of Hamiltonian 
systems and in the semiclassical limit of the nonlinear Schrodinger equation. 

Since S is a hyperplane, T = T a S is independent of a G S. Let W = T x , so any 
(fn, (fin, A) G BV*(N, s, q m , mZ) admits the unique decomposition (f n , (fi n , A) = a + w with 
a G S, w G W, and any U := (u n ,v n ,a) G E decomposes uniquely as U = t + w with 
t G T, w G W. We denote the orthogonal projection maps P : E — >• T, P x : — > 

so that PC/ = t, P ± U = w when U = t + w. We exploit the Hilbert space setting 
and interpret the first variation Vf2f P (fn, (fin, A) as an element of E itself, and project the 
equation Vf2f P (f n , (fin, A) = into the two linear subspaces T and W, 



Pi(r, cr,w) 
P 2 (r, <T,ty) 



P 
P 



BP, 



(7 + W, 



0; 



0. 



(36) 
(37) 



The second equation can be solved uniquely for w = w(r,a) in a neighborhood of S for r 
small, using the Implicit Function Theorem. Because our functional Qf p is smooth we can 
expand w (r, a) in powers of r, and since f2^ p (<7 + u>) is periodic in a the expansion is uniform 
cr. The result below is based on Lemma 2 of ||AmBa 9"8|j , and follows directly from Lemma 



Lemma 5.1 Assume q = q n 



Hp' K 



mZ for some m G N. 



There exist constants r > and S > 0, depending on N, m, k, H, s and a smooth function 



w 



w(r, cr) : (— r , r ) x S ^ W C E 



such that: 



(i) There exists smooth functions W\,W2 such that 

w(r, a) = rwi(a) + r 2 W2(r, a) 
for all \r\ < r and for all a G S; 
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(it) P x [Vfif p {o + w(r, a))] = 0. 

(Hi) Conversely, if P^^Vt^ p (a+w)] = Oforsomer G (— r ,r ) andw G wit/i < 5, 

then w = w(r, a). 

(iv) For any choice of mo, Kq, Sq, H > the constant r may be chosen uniformly for all 
N > 1, 1 < m < mo, 1 < k < k,q, H > Hq and \s\ < sq. 

Parts (i)-(iii) follow easily from the Implicit Function Theorem. The dependences on the 
various parameters in (iv) is more delicate, and was the subject of section 5 of [ AlBeBr PC ]. 



In that section we proved a lower bound on r via a priori estimates on the solutions of (|37|) . 
We note that the interval of validity of the expansion in ||AlBcBr 00|1 was strongly affected by 



the sample width L, with improved convergence with smaller L. In the periodic problem the 
same method as in section 5 of fAlBeBr "00| can be used to obtain the same lower bound on 



the radius of validity, but with the period q = q m replacing L in the expression of the lower 
bound. Since 2q m decreases with increasing field strength H we can expect our solutions to 
have a large range of validity in high fields. 

We define 

S r := {a + w(r, a) : a G S}. 

S r is a smooth manifold parametrized by the hyperplane S. The important role played by 
S r is that it is a natural constraint for Q pp (see Lemma 4 of [|AmBa 9§I1 ), and hence the 
equation (|36|) may be solved variationally: 

Lemma 5.2 Assume q = q m = k = mZ for some m G N. 

(a) If (f n , 0n, A) G S r satisfies D(n PP lSr )(f n , <j> n , A) = 0, then Vn pp (f n , n , A) = in E. 

(b) There exists a constant f = r (iV, m, k, s, H) with < f < r such that if (/„, <fi n , A) G 
BV*(N, s, q m , mZ) is a critical point of Q pp with 

Vt pp {f n , (j) n , A) < 2q m Npr = 2n ™ N r an d | r | < f Q) 

ri 

then (f n ,<fi n ,A) G S r . 

The consequence of this lemma is very important: for r sufficiently small the absolute min- 
imizer will be found by minimization on the finite dimensional manifold S r . We note that 



the parameter dependences for the interval of validity ro in Lemma |5J] do not carry through 
to Lemma |5^. This is because we have no estimate on the neighborhood 5 of Lemma |5j] in 
terms of other parameters, and so we cannot conclude that configurations with very small 
energy must be "close enough" to S to be in the range of Lemma |5.1| (iii). 



The proof of Lemma |5]2| is identical to Lemma 3.2 in [[AlBeBr 00|| and is omitted. 
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6 Vortex lattice solutions. 



We are now ready to treat the finite dimensional problem (|3q) , via minimization of the 
constrained functional Q pp \s r . As opposed to the finite- width case studied in [|AlBcBr 00 1 
the functional Q pp \s r will degenerate at order r, and therefore it will be necessary to carry 
out the expansion to higher order. To accomplish this we must calculate an expansion of 
the solutions of the regular projected equation (|37|) to order r. Fortunately, the regularity 
of Q pp and the Implicit Function Theorem allow us to do the computation explicitly, and 
we will obtain a straight-forward expansion of fl pp \s r which permits direct minimization. 
We summarize our conclusion in the following Theorem: 

Theorem 6.1 Let q m = j^, k = mZ, m G N. 

(i) For every s G R and m = 1,2,..., there exists ro = ro(N, s, K,m, H) > such that 
for all r G (0, ro) the minimizer of Q pp in BV*(N, s, q m , mZ) is a HN^ Sm -periodic solution 
given asymptotically by 



fn = 1 + r 



"2 i " 2(H 2 P 2 +2k 2 ) 



cos(<5 n + Hpx) + cos(5 n+ i + Hpx)) + o(r) 



2 

h{x, z) = H — rjTj cos(<5 n + Hpx) + o(r) 



3t\x) 



2Hp 



&\x) 



(cos(5 n+ i + Hpx) — cos(5 n + Hpx)) + o(r) 
| r K 2 psin(5 n + Hpx) + o(r), 



(3* 



where (S2, 
F{N,s) 



, Sn) is a minimizer of the finite dimensional problem 
inf I cos (^n - <W0 : (£2, • • • , o~n) G R^ -1 , 8x = 0, 5 N+1 = -Hps 



Moreover, the minimum energy is given by: 



inf n pp 

BT»(N,s,q m ,mZ) 



2m qi Np (r + r 2 (C + dF(N, s))) + 0(r 



(39) 
(40) 



where q\ = ■jj- and Co G R ; C\ > are constants independent of N , s,q,m. 

(ii) inf{F(iV, s) : s G R, N = 1,2, . . .} = — 1, and the minimum is attained at (N, s) if 

and only if either: 

N is even and s = 2£qi for £ G Z; 

or if 

N is odd and s = (2£ + 1)^ for £ G Z. 

In either case, the minimizer of Q pp in BV*{N,s,q m ,mZ) is unique and coincides with the 
period-2p in z, period 2q\ in x lattice which minimizes fl pp in BV*{1, q\, qi, Z), satisfying 
( P^ j and given asymptotically by $2$). 
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We observe that the independence of the constants Cq,C\ with respect to N,s,q,m 
means that the dependence of the geometry of II on the energy per cross-sectional unit area 
is entirely encoded in the function F(N, s). In the generic case (i), F(N, s) could have several 
absolute minimizers (62, ■ ■ ■ , Sn), and choosing different sequences of r — > could lead to 
minimizers of Q pp \g r which are determined by different minimizers of F(N, s). In the special 
case (ii) the absolute minimizer of F(N, s) is unique (and non-degenerate) and we obtain 
the smallest possible energy for the Lawrence-Doniach energy, corresponding to the second 
alternative in Theorem 14.11 



Proof of Theorem 6.1: By Lemma 5.1 , we may decompose an element of BV*(N } s, q m , mZ) 
as (f n ,<j> n ,A) = a + w{r,a), with o G S and w G W = [TS] 1 - C E such that the Euler- 
Lagrange equations hold when projected into the space W . Furthermore, we may write 
w(r, a) = rwi(a) + r 2 w 2 (c r , r), and in (/„, <f> n , A x , A z ) coordinates, we denote 

w\ = (un,i,v n ,i,a x ,i,a z>1 ). 



In other words, recalling 

/„ = l+ ru n:1 + 0(r 2 ), 
A x = Hz + ra xA + 0(r 2 ), 



4> n = a n + nHpx + f n ,i + 0(r 2 ), 
A z = ra Zil + 0(r 2 ). 



(41) 



Note also that w\(a) = d r w(0,a). 
Step 1: Expansion of the energy. 

We recognize that the energy can be written in two parts, VL PP {U) = Q PP (U) + rG(U). 
Since fl pp \s r = fl pp (a + w(r,a)) is a smooth function of r and o G S, it admits a Taylor 
expansion at r = of the form, 

2 

nf p (a + w{r, a)) = fi (0) + rfi (1) + -tl^ + 0(r 3 ), 



with 0(r 3 ) remainder uniform in a G S, where: 



0. 



n? p (a + w(r,a)) 



tt {2) = — 
dr 2 



n pp (a) 
d 

G(a) + Wn pp (a)[ Wl (a)}=G(a) 

N \ rx n +2q m 

P/2 (1 — cos(5 n + Hpx)) dx 

n=l Jxn 

2Npq m , 

n pp (a + w(r,a)) 

r=0 
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2WG(a)[w 1 (a)} + V 2 ftf»K Wl ] + Vftf>) 
2VG(a)H + V 2 n^ p (a)[w 1)Wl ] 



d 2 w 
dr 2 



n=l 


+2g m | 


2<i + - 


iV 

+pE / 

n=l Ja; 


v n -\-2q rn 
n 






da x>x 




2 


dz 


dx 





2u n,i + — + - Ox,i(a:, z*)) \ dx 



{(u nj i + Un-i,i)[l - cos(5 n + Hpx)} + sin(<5 n + Hpx)tp njn - lt i(x)} dx 
dz dx. 



Here we denote 



y?n,n-l,l — l>n,l — 



(42) 



2n-l 



Note that the term of order r degenerates- unlike the boundary-value problem treated in 
AlBcBr 00|j the constant phase differences 5 n are not determined at this point, but only at 
the order r 2 ! 

Step 2: Expansion of the solutions. 

To evaluate the next order term we require the first-order correction to the solution 
Wi(a). Implicit differentiation of the equation (j3~7l) with respect to r yields: 



VG(ct) + V 2 n^(a) Wl 



0. 



(43) 



Since V 2 f2^ p (cr) is invertible on W fl43|) determines = Wi(cr) uniquely. 

Using the expansion (|4TD we write the projected equation fl43|) in terms of the above 
coordinates: 



— -u^! + 2u nil = ]- (cos$ nin _i i0 + cos$ n+ i jni0 - 2) , 

= ~[sin(5 n + Hpx) - sm(5 n+ i + Hpx)], 



^dx~ [Vn ' 1 ~ axAX ' Zn 



and 



where 



di{X, z) 



dz" dx 



(44) 
(45) 

(46) 



bi_(x, z) = bf\x), < z < z n , n = 1, . . . , N, 

bf\x) - 6S" +1) (x) = p(v' nl - a Xtl (x,z n )),n= 1,...,N-1, 



d x bi = \ pn 2 sin(5 n + Hpx), 



bf\±L)=Q. 



(47) 
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This system also coincides with the Euler-Lagrange equations for minimizing above in 
the space W. An equivalent way to arrive at these equations is to begin with the equations 



for w(r,a) satisfying (p7\) as derived in section 2.5 of [[AlBeBr 00|| , then take the order r 
terms appearing in each equation. 

These equations may be integrated explicitly: from the current conservation equation 
(|45| ) we obtain: 

k 2 k 2 
v' nl - a Xtl (x,z n ) = C n - -—cos(5 n + Hpx) + ——cos(5 n+1 + Hpx), n = l,...N. (48) 

2Hp 2Hp 

where C n are as-yet undetermined constants. Next we use the equations (47) for the magnetic 
field inside each gap, bi(x) = bi(x, z) for z n _i < z < z n , 

I ^ K 2 

b i ( x ) = -^777 co < 5 n + Hpx) + D n , n=l,...,N, (49) 
Zn 

with D n another set of undetermined constants. Using the jump condition of (|47| ) for the 
magnetic field across each superconducting plane and the periodicity of D n , we see that 

-pC n = (L>„ +1 -AO, n = l,...,N. (50) 

We now solve for the order r term in the gauge-invariant phase y n ,n-i,i (defined in fl42|)) 
using formula (16). We find for n — 1, . . . , N, 

-^-Vn.n-l,l( X ) = ( V 'nl- a xA X > Z n)) ~ Wn-ll~ a xA X > Z n-l)) + Pb^^x) 

ax 

= pD n + C n — C n -\ 

K 2 

+ 2Hp ( cos ^™ +1 + Hpx ^ ~ ( 2 + P 2 ) cos (^" + Hpx) + cos(5 n _i + Hpx)) . 

Since each </? n)n _x,i(^) must be 2g-periodic in x we must satisfy the integrability conditions 
pD n + C n — C n „i = 0. Together with (|50"D, we derive a second order difference equation for 
D n , 

D n+1 - 2D n + D n ^ = p 2 D ni n = l,...,N, (51) 

and note that the e 2 -periodicity condition ([|) of h together with the corresponding condition 
(|35p for 5 n implies D n = D n+ ^. The maximum principle for second order difference equations 
then ensures that D Tl = 0, and hence C n = by (|50"D, and all constants are uniquely 
determined. We may then integrate to obtain: 

k 2 

Vn,n-l,l{ X ) = 2H 2 p 2 ( Sin ^ n+1 + HpX ) _ ( 2 + Sill (^ + H P X ) + sin (^n-l + Hpx)j . 

We note that the arbitrary constant of integration which should normally come with ip n ,n~i,i 
is zero here, since we are solving for each n in the orthogonal to TS. (That constant would 
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result from the order-r correction to our choice of a G S in the minimization problem on 

Sr.) 

In conclusion, the gauge- invariant quantities associated with points (a + w(r, a)) on the 
natural constraint S r are: 

f n = 1 + r -| + 2(g 2p2 2 +2K 2) (cos(5 n + + cos(5„ + i + Hpx)) + 0(r 2 ), 



h(x, z) = H - r-^jT cos(<5„ + Hpx) + 0(r 2 



= r^- (cos(<5 n+1 + Hpx) - cos(5 n + Hpx)) + 0(r 2 ) 



2Hp 



(52) 



j( n \x) = |re 2 psin(5„ + i/px) + 0(r 2 ), 

where (as usual) a = (7(62, ■ ■ ■ , 5n). We emphasize that Lemma [5.1| (i) ensures that all 
remainder terms are uniform in a G S. 



Step 3: Expansion of Q£ \s r - 

We now resolve the degeneracy at order r 2 to determine which choice of o~ = <j(52, ■ ■ ■ , 8 



N 



in S gives rise to stationary solutions of the Lawrence-Doniach system. Substituting and 
computing the integrals, 



Npq r , 



+ Qn 



K 



2H 2 p 2 

K 2 

2H 2 p 2 



2^2 



3 2 _ H 2 p 
2 P H 2 p 2 + 2k 2 



Npq n 



H 2 p 2 



N 



H 2 p 2 + 2k 2 



pJ2 COS ( S n ~ <W), 



n=l 



where we recall q m = mir/Hp. 

In conclusion, we obtain the following expansion of Q pp \s r , with S r parametrized by 
a = a(8 2l ■ ■ - ,S N ): 



Q pp (a + w(r, a)) = 2Npq m jr + r 2 (c o + d ^ cos(5 n - 8 n+1 )) j + 0(r 



(53) 



where Cq G R, C\ > are constants independent of iV, s, g, m. The periodicity conditions 
in BV*(N, s, q m , mZ) carry over to an inhomogeneous boundary condition on S n (see (|35l)), 



^ = 0, 5^+1 = -ifps (mod 2tt). (54) 

We recall once again that the remainder term is uniform in a = (7(82, ■ ■ ■ , 8jy) G 5. 

Step 4: Verifying equation (|40|) , and the conclusion of (i). 

First we observe that by Lemma |4.2| and Lemma |5.2| (b) for < r < f the infimum of 
fl pp in BV*(N, s,q m , mZ) will be attained on <S r . Since <S r is finite dimensional and Q p p \s r 
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is periodic in the local coordinates (82, ■ ■ ■ ,8n), for every < r < r there exists (at least 
one) minimizer, 

ft? p (a r + w(r, <r r )) = inf Q pp \ Sr , a r = a(8 2 (r), 5 N (r)) e S. 
By the expansion (|53|) we have 

Q pp (a r + w(r,a r )) = 2Npq m jr + r 2 (c + d 1 fj cos(5 n (r) - 5 n+ iM)) | + 0(r 3 ) 

> 2iV Mm {r + r 2 (C + C 1 F(iV, S ))} + O(r 3 ), (55) 

since -F(iV, s) is the infimum of the sum of cosines over all possible configurations. To obtain 
a complementary inequality, let (8%, . . . ,8^) be any minimizer of F(N, s), that is 

1 N 



71=1 



(Under the hypotheses of (i) there could be many such minimizer s.) Then, applying ( |53| ) to 
this configuration, we obtain: 

inffif p k < ^ r BP (a* + ™(r,0) 

< 2Npq m |r + r 2 ( C + d i £ ^s(8* n - } + 0(r 3 ) 

= 2Npq m {r + r 2 (C + C 1 F(N,s))} + O(r 3 ). (56) 

Putting together (]55|) and (|56"D we deduce the energy expansion (|40|) stated in Theorem ^TI] (i). 
Moreover, fl56|) implies 



1 ^ 

F(iV, s) < - 2 cos(5 n (r) - 5 n+1 (r)) < F(iV, s) + O(r). 

iV n=l 

Therefore for any sequence of r — > 0, the minimizers ay of fi^ p |5 r form a minimizing sequence 
for the variational problem F(N,s). Hence the a r accumulate as r — > at minimizers of 
_F(iV, s). To be more precise, for any sequence of r — ► there exist subsequences and mini- 
mizers a* = a {5 2, • • • , 8* N ) of F(N, s) such that (along the subsequence) ay — ► cr*. Inserting 
this information into Q52D we obtain Q58p. This completes the proof of part (i) of Theorem 

Step 5: Proof of (ii). 

By the expansion fl40|) the problem reduces to determining for which lattice parameters 
N, s does F(N, s) attain its lower bound of —1. This lower bound is achieved if and only if 
the boundary condition (|5l) admits a choice of 8 n with 8 n+ i — 8 n = ir (mod 2n.) 



S. Alama, A.J. Berlinsky, & L. Bronsard 



31 



When s ^ -^Z the lattice is frustrated since the space BV*(N, s, q m , mZ) does not admit 
the configuration with 8 n+ i — 8 n — ir (mod 2%.) In that case we obtain F(N, s) > —1, and 
the energy per unit area of the minimizer in BV*(N, s, q m , mZ) will be strictly larger than 
this absolute minimum value, for all sufficiently small r > 0. 

If 

N is even and s = 2£^ for £ e Z (57) 

or if 

N is odd and s = (2£ + 1)-^ for £ E Z, (58) 

the choice 8 n = (n — l)n (mod 2tt) is allowed by (0) and the infimum F(N,s) = —1 is 
attained. Define j:Rx R^ -1 — > R with 

9(r5 2 Sn) = 1 f ^ P ^ + w(r,a)) \ 
'' r 2 \ 2mqiNp J 

1 N 

= Co + d— J2cos{8 n -8 n+1 ) + 0(r), 

iV n=l 

by (|53|) . In particular, g is smooth and g(0, 82 ■ ■ ■ , 5jy) is minimized if and only if 8 n = 
5* = (n — l)n mod 2n. Now, (82, • • • , 8%) is a non- degenerate minimizer of g(0, 82 ... , 5at): 
its Hessian is the familiar tridiagonal matrix with 2 on the diagonal and —1 on each off- 
diagonal, associated with the (positive-definite) second-order difference operator. By the 
Implicit Function Theorem we conclude that for all sufficiently small r > the function 
g(r, 82 ... , 8n) has a unique minimum at (^(r), . . . , 8jsr(r)), with 8 n (r) = 5* + 0(r). In other 
words, setting a r := cr(5 2 (r) . . . , <5jv(r)), 

inf^ p | 5r = af P K. + u7(r,a r )). 



By Lemma 5^ when < r < fo this gives the global minimizer of Q pp in BV*(N, s, q m , mZ), 
and inserting these optimal values for 8 n into the asymptotic formulae (|5^) we obtain (|29|). 
This proves that the minimizers for any period geometry satisfying (|57D or (p8|) coincide up 
to order r. We will do better, and show that they are actually identical. 

To gain a complete understanding of the absolute lowest energy solution we first consider 
the special case N = 1, s = q% = 71 /Hp, and m = 1. In this case, there is no degenerate 
manifold: the r = problem has a unique, non-degenerate solution and the perturbation 
is regular. In particular the second-order term in the energy expansion (^) is completely 
determined by (|54f) : since 8 = and (for N = 1) 8n+i = 8 2 = —Hsp = —n we have 
F(l, q\) = — 1. We obtain for all sufficiently small r > a unique solution 0* , A 1 ) which 
minimizes fl pp in BV*(1, s, qi, Z). The gauge-invariant quantities associated to (/*, 0* , A 1 ) 
will be 2gi-periodic in x and shifting z by p results in a translation by a half-period q± in x 
(see (HD.) 
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Now return to the cases (|57| ) or where N > 2 and F(N,s) = — 1. We denote by 
<5* = (n — l)7r mod 27r, the unique absolute minimizer of F(N, s), and by (S 2 (r), . . . , <5jv(r)) 
the coordinates of the absolute minimizer (f n , <fi n , A) G BV*(N, s, q m , mZ) of fif p , which we 
know lies on S r for all small r, and for which 5 n (r) = <5* + 0(r). 

It is easy to verify that the £vP*(l, gi, q±, Z)-minimizer (/^, 0* , A 1 ) is also ILv^^-periodic 
when iV, s satisfy (f>% or fl58|) , and in fact (f^^A 1 ) solves the Euler-Lagrange equations 
for Q pp in BV*(N, s, q m , mZ). Moreover, (fn,4>n, A 1 ) also satisfies the expansion of energy 
given by ((EJ) on the period parallelogram rLv,s l(?m , and therefore by Lemma |572| it lies on the 
constraint manifold S r C BV*(N, s, q m , mZ) in cases fl57|) and (|58|) . Finally, the coordinates 
of 4^, A 1 ) on S r C BV*(N, s, g m , mZ), which we denote by (5^(r), . . . , 5^\r)), satisfy 

5 { n ] =5: + 0(l). 

Since (#3, . . . , 5jy) is a non-degenerate minimizer of F(N,s) we must have <5 n (r) = S^(r) 
for all sufficiently small r, that is (/„, <f) n , A) = 0*, ^l 1 ) for all small r, that is the 



minimizer s in BV*(N, s, q m , mZ) with 
^(1,^,^,2). 



or ( |58D coincide exactly with the minimizers in 



This concludes the proof of Theorem 6.1 





Finally we prove Theorem O] , In case g ^ J^Z or fc„ ^ mn (for m G N constant) we 
have by Lemma fO 



e r (iV, s, q,k) := inf Q 

BV,(N,s,q,k) 



BP 



> inf 

BV t {N,s,q,k) 



n BP 



=: u; (iV, s, g, fc) > 0, 



with uq(N, s, q, k) a constant independent of r. Using the expansion (^) of the minimizing 
solution in BV*(l, jj-, jj-, Z) we obtain the first alternative in the statement of the Theorem. 



When g 6 -^Z the remaining statements follow as a corollary to Theorem |6TT . 





Remark 6.2 The result contained in Theorems 3~T and O, namely that the smallest possi- 
ble energy per unit cross-sectional area is obtained with N = 1, s = qi, m = 1, confirms the 
prediction of a period- 2p in z staggered lattice solution made by Bulaevskh-Clem [[BuCm 91 



In that paper the authors indeed claim that it minimizes energy among competing configu- 
rations, but only one other vortex lattice is treated in their paper (a period-4 lattice), and 
no indication is provided as to how they deduced the geometry of their solution. Here we 
have shown much more: we know that it is the minimizer among all periodic solutions, in 
the regime r << 1. 

The "vortex plane solution" of Theorodakis-Kuplevaksky ( ||Th 90|| , |[K 99|| ) is obtained by 
taking 5 n = for all n (when the value of s permits such a choice.) Note that such a choice 
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maximizes on the constraint set <S r . Since these would also constitute non-degenerate 
critical points of the finite-dimensional function g(r, 8 2 , ■ ■ ■ , 8n), by the above arguments 
they describe bona fida solutions to the Lawrence-Doniach system with periodic boundary 
conditions, but they are unstable. 

Remark 6.3 When s 7^ jq±, j G Z the trivial choice of minima 5 n — <5 n _i = ±7r is not 
admissible, and the lattice is "frustrated". As previously mentioned, in this generic case 
the minimizer of F(N, s) might be non- unique, leading to different asymptotics for the 
minimizing solutions along different subsequences r — > 0. Nevertheless, the dependence of the 
energy minimizers on subsequential limits r — > could be eliminated when it is known that 
the absolute minimizer (8 2 , . . . , 5^) of F(N, s) is unique and non-degenerate. An example is 
when the minimizer of the finite dimensional problem ( |39"D (8%, . . . , 8^) satisfies 

C n := cos(<5; - 5* n+1 ) < for all n = 1, . . . , N. (59) 

This will be the case when the parallelogram TlN,s,q is very close to the optimal ones described 
by fl57p and flSBp. Assuming ( |5^) holds for the minimizer, a simple calculation shows that 
the Hessian is the ((JV — 1) x (N — 1)) tridiagonal, symmetric matrix D 2 g(0, 8%, . . . , 8^) =: 

[■Wm,n]m,n=2,— ,iV with 

M ntn+1 = C n n = 2,...,JV-l; M n , n = -(C n _i + C n ), n = 2,...,N. 
A null vector v — (v 2 , ■ ■ ■ , v n) of M satisfies: 

C n {v n+ i - v n ) - C n ^i{v n - = 0, n = 3, . . . , N - 1; 

C2<>3 - f 2 ) - CiW 2 = 0, C N V N - C N ^{V N - V jv-i) = 0. 

If i>2 = then clearly u = 0, so we may assume v-i > 0, in which case equations n = 
2, . . . , N — 1 imply < v 2 < . . . < v^. But this contradicts the n = N equation, and 
therefore the only solution is the trivial one. In conclusion the minimizer is non-degenerate, 
and we can repeat the same arguments as in the cases ([57]) and fl5"8D to conclude uniqueness 
for energy minimizers. 

7 The periodic finite layer case 

Finally, we consider the case of a finite number of planes N, each of infinite extent in x and 
y, assuming that the currents and field strength are periodic functions in x. Since this case is 
very similar to the doubly periodic case treated in the previous sections we give an outline of 
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f q 

/ (pi(x)dx; 

J-n 



how to modify the formulation of the problem and its solution to fit this somewhat simpler 
case. 

In the finite-layer case the 't Hooft condition is greatly simplified, since (by the argument 
of Theorem j37T|) we may take A to be periodic in x. Let k = (k n ) ne z with k = 0, k ^ 0. 
We say that (/„, <p n , A) belongs to the periodic class V = V(q, k) if there exists a constant 
lj G R such that: 

A6C(R 2 ;R 2 ), U e H l per (R), <f> n e #/ 0C (R); 

<f> n (x + 2q)=(f> n (x)+u + 2Trk n , n = 0,...,JV; (60) 
ri ri 

-g ->-i 
A(x, z) = (Hz, 0) + -Q, with e G ^ 2 oc (R 2 ) 

e(x + 2g,2 ; )=e(x,2:), 0) = ^(x, JVp) = 0. , 

If A; = we define V(q,0) by omitting (|6T|). Following the proof of Theorem |3.1| , any 
configuration satisfying a single 't Hooft condition (in x with u x ) as in ([5]), @ is gauge- 
equivalent to an element of V(q, k). Furthermore, in this class the L 2 norm of curl A controls 
Am H 1 . 

As in the doubly periodic case, the observables are determined entirely by the values of 
{fni (fin, A) in a single period x G [—q, q), but now the number of planes is finite (indexed by 
n = 0, . . . , N) and each is independent. Hence we define the Gibbs free energy by integration 
over a single period and summation over the iV + 1 planes, 

N 



(61) 
(62) 



J - ( i n=0 



« n -A x {x,z n )ffi 



dx 



r ri N / \ 

+ 2 j_ P £ {f" + fn-l ~ Vnfn-l COS^^)) dx 



-1 n=l 

I r N p fi (dA x 

q 



K JO 



dz 



dAz 
dx 



— H ] dx dz. 



By taking variation of Q^ er (f n , <p n , A) in the subclass of periodic functions, we see that the 
Euler-Lagrange equations are exactly the same as for the finite width sample case, see (2)- 
(5), (7) of [|AlBeBr 00|| . That is, they coincide with (|Il])-(|15|) for n = 1, . . . , iV — 1, and the 
equations involving the top and bottom surfaces n = 0,N must be modified to reflect the 
fact that these two planes have only one "nearest neighbor" . 



We may now continue as in the doubly periodic case. By Lemma |4.2| the minimum 
energy will be of order r if and only if we choose q = rrij^ and k = mZ, with m = 1, 2, . . .. 
Furthermore, with these choices the minimizers will have f n = \ip n \ ~ 1, and the use of polar 
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coordinated for the order parameter is well justified. We choose q = k = Z, and write 
p # : = V(j^, Z) in the following. 

For the number N of planes, k, and H fixed, define the free energy per unit area in a 
period strip of width 2q with winding numbers k by: 

e r (q,k) := inf {^(/n, <!>n,A) : e%,fc)}. 

We obtain the following result: 
Theorem 7.1 Let iV, iJ, k be fixed. 

(a) If there exists m e N sitc/j £/icrf 

Q = Qm '■= ~rr~ an " ^ = (63) 
Hp 

then there exists r = r (N, H, K,m) > swc/i that for all < r < r , e r (q,k) = 
e r (-j^,Z), and the minimizers of Q^ er in P(g m , mZ) coincide with the minimizers of 
ni er inV*. 

(b) For any other choice of q, k there exist constants r 1; oj\ > (depending on N, H, K,q,k) 
such that e r (q, k) > uj\ for all < r < r%. 

We note that, as in the bi-periodic case, we may obtain an expansion of the minimizing 
solution in powers of r. The minimizer in the special space V* (which gives the absolute 
minimum of energy per unit area) coincides with the period- 2p in z lattice found in the 
bi-periodic case at order r, except for an edge effect in the order parameter in the top and 
bottom planes. More precisely, the fields and currents h(x, z), j^ n \ j( n ' in a finite stack still 
satisfy (p9|) for each n = 0, . . . , N, and /„ coincides with the expression given in ( p9l) for 
n — 1, . . . , N — 1 but is modified (see (^) and (|69D below) at order r for n = and n = N . 



The proof of Theorem [fj] follows almost line-for-line the degenerate perturbation proce- 
dure of the previous sections. In particular, the minimum value of Q^. er at r = is attained 
by elements in the (N — l)-dimensional hyperplane 

S := {(/„, <f) n , A) EV* : f n = 1, <f) n (x) = a n + nHpx, (64) 
A x = Hz, A z = 0, where «o = ai = 0, and «2, . . . , oin G R.} 



The only significant difference with the doubly periodic problem occurs at the top and 
bottom layers, n = 0, N. The first-order expansion of the solution in r requires that these 
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be treated differently. For example the equation of current conservation ( pE5|) now gives 

{2 2 
C n ~ 2% cos(5 n + Hpx) + ^ cos(5 n+ i + Hpx), n = 1, . . . N - 1; 
^0+ 2^ c os(^i + Hpx), n = 0; 

C n ~ ^cos(5 N + Hpx), n = N; 

(65) 

The order r term in magnetic field &i is exactly as in the doubly periodic case and is given 
by (P&D (with as yet undetermined constants D\, . . . , -Djv-) Therefore, using equation (ffl) 
and the jump conditions (p7[), we obtain the following conditions on the constants C n and 
A,: 

A t+ i - 2D n + L> n -i = P 2 A, n = 1, ■ ■ ■ , N - 1, (66) 
D 2 -2D 1 =p 2 D 1 , D N _ 1 -2D N =p 2 D N} (67) 
-pC n = {D n+1 -D n ), n = 1, . . . , N - 1, 
-pC = A, VC N = D N . 

The maximum principle for the second-order difference equation (|66| ) with boundary condi- 
tion fl6T|) implies that the unique solution is D n = 0, C n = for all n. In consequence the 



perturbed manifold S r consists of the same configurations (|52 ) as for the bi-periodic case, 
except for the superconducting order parameter which coincides with the expression given 
flS"^ ) for n = 1, . . . , N — 1, but for n = or n = iV we obtain 

= 1 + \ (4 + 2(flv+2^) C08(fe + + ° (r2) - (68) 

A = 1 + 5 ("5 + ^ cos(i/pl) ) + ° (r2) - (69) 

Substitution of fl5"2|), (|68|), (|69|) into f2^ er leads to an expansion of the energy in the same 
form as (|53), 

Qf> + w(r, a)) = 2Npq jr + r 2 ^ + C 1 1 £ cos(5 n - <Wi)l } + 0(r 3 ), 

except for the constants Co G R and C\ > which differ from the doubly periodic case due 
to the slightly different form of solutions for the top and bottom planes. The significant 
difference from the previous case is that when there are finitely many planes there is only 
the single constraint 5± = (which comes from removing the translation invariance in the 
definition of the space V*). By the same arguments as in the previous section the minimizer 
of QP er will be determined by minimizing the leading term in the energy expansion, 

N-l 

G(5 2 , ...,S N )= cos ($n - <Wi), 5i = 0. 

n=l 
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By inspection, the minimizer is obtained by choosing 5 n+ i — 5 n = ±7r; for example 8\ = it, 
8 2 = 0, 83 = 7i, . . . Unlike the doubly periodic case we are always free to make this choice, 
and it is easy to verify that this configuration gives a non-degenerate minimizer of the finite- 
dimensional function g(0, a), a = (62, ■ ■ ■ , 8n)- In particular the period- 2p in z, period-2gi in 
x lattice appears naturally (at order r) in the solutions. (See figure 1.) As in the finite-width 
case the solutions feel the top and bottom edges only in the first plane and first gap at order 
r; expansion to higher orders will reveal the effect of the finite stack at order r k in the kth 
stack from the top or bottom. 



8 Conclusions 

Finally, we summarize our results on the periodic problem, and compare with the (quite 
different) conclusions obtained for finite- width samples in our previous paper ||AlBeBr 00|| . 

As we have seen, in the limit r — > energy minimization selects a preferred period 
geometry and quantized flux from all possible periodic configurations represented by the 
spaces £vP*(iV, s,q,k). The optimal solution is -^-periodic in x, and repeats itself with a 
horizontal shift of a half-period in x when we climb from one superconducting plane to the 
next. Each period parallelogram contains exactly one quantum 2tt of flux, and one Josephson 
vortex per period in each gap between the superconducting planes. The Josephson lattice 
geometry is the same as was predicted by Bulaevskii & Clem |BuCm 91|] , but with our 



rigorous analytical approach we may now assert that it is the unique energy minimizing 
periodic configuration (among all possible geometries) for all sufficiently small r. 

In the finite-width case [|AlBcBr 00|| the conclusions were surprisingly very different. For 



any finite width sample, — L < x < L, when sin(HpL) 7^ the unique energy minimizer for 
r ~ is a vortex plane configuration, with Josephson vortices vertically aligned, and the 
magnetic field approximately uniform in z (except for edge effects at the top and bottom of 
a sample of finitely many planes.) The exceptional values of the applied field H for which 
sm(HpL) = correspond to first-order phase transitions occuring when a new vortex plane 
is nucleated into the sample from the lateral edges. 

We can explain this apparent conflict by examination of the expansion of the energy 
in powers of r near the degenerate manifold S. In (23) of [ [AlBeBr 00| | vortex planes are 



preferred at order r in the energy because of a surface term, a quantity which scales like the 
length Np of the lateral edges of the sample cross-section. In the periodic case this term 
does not appear, and in (^) the distinction between lattice geometries appears at order r 2 , 
in a term which scales as the cross-sectional area 2qNp of the sample. For any size sample, 
by making r small enough the order r surface term will dominate, but the value of r must 
decrease if the surface term is to continue to prevail with increasing sample width L. If we 
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Figure 1: Period-2p (in z) vortex lattice, for a sample with a finite number of superconducting 
planes (Indicated by horizontal dotted lines.) Horizontal arrows indicate the in-plane currents 
and the vertical arrows depict the Josephson currents ji n ^ between adjacent planes. The magnetic 
held h(x,z) and supercurrents j x , j z are periodic with period 2p in the z-direction (and period 
in x). The vortices (local maxima of h) lie along the starred segments, and form a staggered 
lattice. If we choose the midpoint of each segment to label each vortex the resulting lattice is 
diamond-shaped. 

try to keep r fixed while increasing L then inevitably the order r 2 term will compete with 
the order r term. At that point the perturbation expansion will surely have lost its validity. 

This general argument is supported by the analysis of the range of validity of the expan- 
sion in Lemma |5.1|, presented in section 5 of ||AlBeBr 00|| . Indeed, it seems clear that the 
radius of validity of the expansions in r deteriorates with increasing sample width L. (See 



Remark 5.3 of ||AlBeBr 00|| .) 

If we approximate a macroscopic sample by an infinite one and seek periodic solutions, 
the period plays the role of the sample width in the estimates of the interval of validity of 



section 5 in [[AlBeBr 00|| . Since the period of the absolute minimizer is given by we 



can expect the interval of validity to extend to physically appropriate values of r when the 
applied field H is large enough. Therefore we may apply our analysis to the transparent 
state of the high-T c superconductors in high external fields. 
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